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Abstract 

Recent  technological  advancements  in  solar  thermal  rocket  propulsion  and  solar  orbit 
transfer  vehicles  make  it  critical  to  perform  additional  engine  performance  analyses.  Several 
system  level  flight  demonstrations  are  imminent.  Space  flight  hardware  component  testing  is 
being  conducted  at  the  Air  Force  Research  Laboratory,  Edwards  AFB,  California.  The  focus  of 
current  research  is  engine  and  nozzle  configurations  for  a  solar  orbit  transfer  vehicle.  The 
optimal  design  must  produce  1-10  pounds  thrust,  perform  at  high  Isp  and  be  compatible  in  a 
hybrid  of  spiral,  perigee,  and  apogee  (multi-burn)  configurations.  The  nozzle  material  must  not 
ablate  when  subjected  to  extreme  thermal  loading,  yet  be  durable  enough  to  withstand  widely 
varying  temperature  differentials  during  frequent  thermal  cycling.  This  paper  addresses 
propulsive  needs  in  the  orbit  transfer  arena  and  defines  governing  upper  stage  vehicle  engine 
equations.  These  equations  are  modified  versions  of  rocket  engine  equations  used  for  chemical 
systems.  The  correction  factors  and  modifications  are  for  Solar  Thermal  Propulsion  specific 
hardware. 


Introduction 

Solar  Thermal  Propulsion  (STP)  uses  large  parabolic  off-axis  mirrors  to  focus  sunlight  into 
an  absorber-heat  exchanger.  The  solar  concentrators  transfer  a  minimum  of  3000  K  into  the 
propellant  (Fig.  1).  A  blackbody  absorber  is  a  heat-transferring  engine  that  will  not  emissively 
degrade  the  incoming  solar  flux.  Gaseous  hydrogen  propellant  is  directed  from  an  inlet  and 
distributed  throughout  the  heat-exchanging  cavity.  Hydrogen  gas,  used  for  its  low  molecular 
weight,  must  be  channeled  through  the  absorber  walls  because  it  is  largely  transparent  to 
sunlight,  and  will  not  absorb  heat  on  its  own.  This  heat  transfer  process  allows  the  hot  hydrogen 
to  be  thermodynamically  expanded  through  a  nozzle  located  aft  of  the  engine  absorber.  This' 
propulsion  concept  is  simple  and  effective. 

The  current  fleet  of  spacecraft  in  the  US  inventory  can  place  less  than  1%  of  their  total 
launch  mass  into  high  orbits.  The  United  States  Air  Force  is  continually  exploring  innovative 
propulsion  concepts,  such  as  STP,  to  place  military  payloads  into  space  reliably  at  lower  cost. 
Ground  demonstrations  have  proven  the  capability  of  engine  performance  and  concentrator/ 
strut  deployment.  An  upcoming  flight  test  experiment  will  validate  STP’s  performance  potential 
as  a  reliable  orbit  transfer  vehicle. 
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Figure  1.  STP  System  Hardware 


A  Low  Earth  Orbit  (LEO)  to  Geosynchronous  Earth  Orbit  (GEO)  transfer  vehicle  (upper 
stage)  must  exceed  current  vehicle  performance  levels  if  it  is  to  replace  traditional  (chemical  and 
electrical)  rocket  systems.  Although  high  specific  impulse  (Isp)  reduces  propellant  mass,  it 
compromises  thrust  levels.  Trade-off  studies  prove  that  the  propellant  storage  volume  saved 
would  enable  more  payload  to  be  manifested  in  a  launch  vehicle  shroud  with  only  a  moderate 
increase  in  trip  time.  For  a  typical  solar  orbit  transfer  trip  time  would  be  30  days  for  a  multi-bum 
transfer,  utilizing  a  hybrid  of  perigee  and  apogee  bums.  To  effectively  reach  the  GEO  altitude, 
the  spacecraft  must  perform  roughly  10  times  as  many  perigee  as  apogee  burns  (Fig.  2). 


Perigee  Burns 
Roughly  90% 

of  Total 


n-^rri  n+ 


Apogee  Burns 
Roughly  1  0% 
of  Total 


Figure  2.  Transfer  Type  for  30-Day  Trip  Time 


The  Basic  Solar  Engine 


The  basic  solar  engine  offers  substantial  performance  efficiency  improvements  over 
traditional  chemical  rocket  performance.  Chemical  rockets,  which  have  the  energy  source 
coincident  with  the  propellant,  are  fundamentally  limited  in  their  achievable  specific  impulse  by 
the  strength  of  the  chemical  bonds  of  the  propellant  (Isp  =  450s).  If  high  specific  impulses  are  to 
be  achieved,  an  energy  source  other  than,  or  in  addition  to,  the  propellant  must  be  used.  The 
solar  concentrators  serve  this  purpose  for  a  STP  without  having  to  use  combustion. 

The  two  most  important  rocket  engine  performance  parameters  are  thrust  and  Isp.  Thrust  is 
the  final  velocity  of  the  propellant  times  the  mass  flow  rate  out  of  the  nozzle.  The  solar  engine 
operates  by  expanding  hot  hydrogen  gas  and  accelerating  it  to  produce  thrust.  Thrust  is  then 
divided  by  the  total  propellant  mass  flow  rate  to  determine  Isp.  Using  basic  rocket  propulsion 
equations  we  can  solve  for  these  parameters: 


P thrust =  1  /2  mdot  Vexit  ("I ) 

Equation  1  shows  how  the  mass  flow  rate  is  related  to  power  (Pthmst  )>  while  Equations  2  and  3 
present  the  basic  definitions  of  Isp  and  thrust. 

F  =  mdot  Vexit  +  m  vdot  exit  or  F  =  mdot  vexit  +  APe*it  Aexit  (2) 

Isp  (g)  =  F  /  mdot  or  Pthrust  =  1/2  F  Isp  (g)  (3) 


These  equations  work  fairly  well  for  defining  ideal  rocket  engine  performance;  unfortunately 
several  key  and  uncommon  factors  must  be  addressed  which  play  a  significant  role  in  designing 
an  optimal  solar  engine.  For  example,  boundary  conditions  and  thermal  losses  for  flowing  hot, 
corrosive  hydrogen  through  a  blackbody  absorber  need  to  be  considered,  and  the  governing 
equations  must  be  modified  to  accommodate  these  changes. 

Why  Basic  Rocket  Equations  Fail  Solar 

As  propellant  hydrogen  is  lightweight,  cryogenically  storable,  and  ideal  for  STP  propulsive 
needs.  However,  when  heated  to  extremely  high  temperature  (3000  K),  dissociation  of  diatomic 
hydrogen  molecules  becomes  a  concern.  From  a  thermodynamic  standpoint,  the  rate  at  which  H2 
is  dissociated  into  H  will  not  affect  the  ratio  of  specific  heat  (y).  This  is  because  the  rate  of 
dissociation  is  sufficiently  small  enough  at  these  temperatures  that  the  ratio  of  specific  heats  is 
essentially  constant. 

When  the  propellant  is  expelled  through  a  rocket  nozzle  with  no  further  chemical  reaction, 
the  flow  is  said  to  be  “frozen”.  If  frozen  flow  is  assumed  (for  traditional  propulsion  equations),  the 
mixture  of  these  two  gases,  diatomic  and  monatomic,  will  not  change  as  it  is  accelerated  through 
the  engine.  However,  more  accurately,  the  energy  used  to  disjoin  the  hydrogen  is  conserved  by 
the  recombination  of  H  into  H2  along  the  nozzle  surfaces.  Although  this  operation  returns  some 
of  the  energy  to  the  gas,  it  creates  a  boundary  layer  thrust  loss  and  reduces  Isp. 

The  blackbody  radiator  employed  by  STP  raises  some  important  performance  analysis  and 
sizing  concerns.  Thermal  losses  and  power  drops  must  be  accounted  for  in  new  thrust  and  Isp 
equations  to  optimize  the  engine  design.  The  correction  factor  that  accounts  for  power  losses, 
xett,  is  a  measure  of  efficiency.  Total  efficiency  is  a  sum  of  the  internal  efficiency,  rijnt,  and  the 
propulsive  efficiency,  qp.  For  solar  propulsion  qjnt  is  assumed  to  be  0.50.  This  factor  can  be 
placed  within  the  governing  rocket  equations.  Also  the  effective  exhaust  velocity,  c,  is  used 
instead  of  the  ideal  exhaust  velocity,  veXjt,  to  incorporate  these  losses.  The  efficiency  of  the 
engine,  teir.  encompasses  thrust  vector  loss  along  the  rocket  axis  (quantified  by  the  cosine  of  the 
nozzle  angle)  and  overall  nozzle  power  loss.  The  nozzle  angle  is  defined  in  Figure  3  as  alpha. 
The  nozzle  used  for  STP  is  a  converging/diverging  nozzle. 


Equations  4-7  account  for  losses  that  are  generally  disregarded  during  the  design  process 
for  chemical  systems.  This  is  not  to  suggest  that  chemical  propulsive  systems  do  not  encounter 
significant  performance  drops  due  to  power  losses;  quite  the  opposite,  chemical  systems  must 
deal  with  possibly  greater  net  losses  due  to  the  complexities  associated  with  combustion.  While 
the  inclusion  of  efficiency  terms  may  lead  to  a  more  closely  traceable  STP  system,  the  values 
attributed  to  these  terms  must  be  analytically  determined  and  then  validated  through 
experimental  testing.  It  is  through  extensive  “ground”  engine  testing  that  actual  vs.  real  thrust/ 
Isp  levels  can  be  quantified. 


Teff-  hint  hp 

(4) 

Pthmst  =  Teff  Ah  mdot/p  =1/2  F  Isp  (g) 

(5) 

1/2  F  Isp  (g)  =  1/2  mdot  c2 

(6) 

Isp  =  1/g  (2  tetr  Ah  mdot  /  p) 

(7) 

The  nozzle  power  loss  is  defined  as  5%.  Therefore,  the  effective  velocity  used  in  Equation 
6  can  be  solved  using  Equation  8.  The  subscripts  used  to  denote  location  along  the  thrust  axis  of 
the  engine  are  shown  in  Figure  3.  The  effective  velocity,  c,  is  defined  by: 

c  =  v2  +  (P2  -  Ps)  A2/mdot.  (8) 


ALPHA  <15deg 
BETA  <  60deg 


Figure  3.  Subscript  Location  and  Characteristic  Angle  Design  Criteria 


Change  in  enthalpy,  Ah,  must  include  thermodynamic  losses.  The  parameter  Ah  =  habsorber  exit  - 
^thruster  exit-  The  heat  of  formation  is  defined  as  the  positive  interaction  required  to  form  the 
original  composition  from  its  products  (H  and  H2)  at  a  constant  pressure  and  prescribed 
temperature.  For  the  STP  case,  the  heat  of  formation  refers  to  the  amount  of  dissociated 
hydrogen  present  in  the  total  propellant  mass  as  a  result  of  expansion  and  heating.  In  Equation  9 
this  is  represented  by: 


Ah  =  Sn2Ah2°  -  Sn3Ah3°  (9) 

where  the  subscripts  2  and  3  describe  whether  this  value  is  taken  prior  or  after  full  expansion 
within  the  engine  (leaving  the  absorber  or  leaving  the  thruster),  n  is  the  number  of  moles  and  Ah 
is  the  heat  of  formation  per  mole. 

Solving  Equation  9  yields  the  following: 

1 . )  Clarifies  why  real  vs.  theoretical  number  divergence  occurs, 

2. )  Enables  more  reliable  thermal  modeling,  and 

3. )  Reduces  the  number  of  iterations  required  to  produce  acceptable  performance 

levels  for  future  designs  that  meet  mission.  These  equations  set  the  baseline  for  STP 
design  criteria. 


The  Design  Process 


This  section  presents  the  detailed  equations  used  to  design  our  STP  engines.  Since  our 
engine  does  not  require  combustion,  only  high  temperature  losses  need  to  be  considered.  There 
are  no  combustion  or  chemical  losses  present  in  a  STP  system.  Since  the  designer  is  aware  of 
the  performance  values  (Isp  and  T)  needed  to  meet  mission  requirements,  designing  a  STP 
engine  is  a  simple,  executable  process. 

To  determine  nozzle  size  used  for  given  engine-operating  conditions,  important 
temperatures  and  pressures  must  first  be  calculated,  and  then  measured  experimentally. 
Equations  10-13  theoretically  determine  the  temperatures  and  pressures  extracted  from 
chamber/absorber  conditions  related  to  the  nozzle  throat.  The  theoretical  values  for  these  actual 
conditions  should  be  confirmed  experimentally.  The  solar  laboratory  facility  supplies  the  means 
to  measure  the  calculated  values  using  thermocouples  and  pressure  transducer  input  to  data 
acquisition  hardware/software.  This  allows  the  experimenter  to  actively  monitor  the  test  engine 
and  insure  that  it  is  operating  under  its  design  parameters. 

T2  is  the  temperature  of  the  hydrogen  gas  at  the  nozzle  throat.  It  is  given  by  Equation  10: 

T2  =  Ti[1/(1+(Y-1)/2)].  (10) 

The  ratio  of  specific  heats,  y,  is  equal  to  the  nondimensional  value  1 .4  (even  when  the  hydrogen 
propellant  is  subjected  to  extremely  high  temperatures).  This  yields  a  throat  temperature  of: 


T2  =  0.833  Ti  (11) 

where  Ti  is  the  absorber  chamber  temperature  given  in  degrees  Rankine  (°R).  With  an  absorber 
temperature  of  3600°R,  the  throat  temperature  then  can  be  derived. 

Pt  is  the  hydrogen  gas  pressure  at  the  throat.  The  pressure  at  this  location  is  less  than  the 
chamber  pressure,  due  to  the  acceleration  of  the  gas  through  the  engine.  The  value  for  Pt  is 
given  by: 


P2  =  Pi  [1+ (y-2)/2]'r/(r1).  (12) 

The  chamber  pressure  is  typically  set  at  150  psia  for  STP  applications.  This  value,  g,  was 
experimentally  validated  in  conditions  that  met  the  limitations  and  mass  flow  requirements 
defined  in  Equation  6.  Inserting  this  value  for  g  into  Equation  8  gives: 

P2  =  0.528  P,.  (13) 

To  obtain  maximum  thrust,  the  hot  hydrogen  gas  must  be  expanded  in  the  diverging  section 
of  the  nozzle.  Since  the  amount  of  energy  needed  to  accelerate  the  hydrogen  creates  a  net 
pressure  drop  through  the  nozzle,  the  next  step  in  the  design  process  is  to  locate  the  exact 
location  where  the  gas  pressure  is  equal  to  the  local  atmospheric  pressure.  This  area  which  is 
denoted  by  the  subscript  3  in  Figure  3  is  commonly  called  the  nozzle  exit  area. 

First,  to  define  the  exit  area  of  the  nozzle,  the  Mach  number  must  be  found.  The  Mach 
number  is  the  ratio  of  the  velocity  of  the  gas  to  the  local  speed  of  sound.  The  perfect  gas 
expansion  expression  can  be  used  since  our  engine  involves  no  combustion. 

M32  =  2/(y-1)t(p1/p3)(r1)/Y-1]  (14) 

P3  iS  the  value  of  the  pressure  at  the  nozzle  exit  plane.  For  a  STP  system  in  LEO,  a  good 
approximation  of  patm  is  0.2  psia  due  to  plume  interaction  or  backflow.  The  nozzle  cannot  be  fully 
expanded  and  its  length  must  be  truncated  to  discourage  flow  separation  along  the  nozzle  walls, 
and  subsequent  shock  development.  It  has  also  been  numerically  determined  that  an  overall 
nozzle  area  ratio  (exit  to  throat)  be  assumed  to  be  50:1.  Once  the  Mach  number  is  calculated, 


determining  the  geometric  qualities  of  the  operating  system  can  begin.  The  exit  area  and  the 
throat  area  are  given  in  Equation  15: 


A3  =  A2/M3  [1  +  ((y-1)  /  2)  M32  /  ((y-1)  /  2)] (rM) ' 2(r1)  (1 5) 


where 


A2  =  mdot  /  p2  (R  T2/  (y  gc))1/2.  (16) 

Solving  for  mdot  using  the  corrected  Equation  7  yields  a  propellant  flow  rate  of  0.002  lb.  per 
second.  The  throat  area  is  then  calculated.  Using  Equation  10  we  find  the  exit  area.  Note  that 
the  area  expansion  ratio  is  usually  prescribed  to  be  50:1.  Using  typical  STP  operating  pressures 
and  temperatures  (150  psia  and  3600°R)  the  actual  expansion  ratio  is  32:1.  Solve  Equations  10 
and  1 1  to  calculate  the  diameter  of  the  orifices.  While  these  equations  take  into  consideration 
power  loses  due  to  propellant  and  propulsive  inefficiencies;  they  do  not  account  for  boundary 
layer  accumulation  or  flow  divergence.  The  next  section  will  deal  with  boundary  layer  effects 
and  characteristic  lengths  of  the  engine.  Since  the  nozzles  are  designed  to  be  underexpanded, 
there  is  no  accumulation  of  shock  waves  within  the  cone  of  the  diverging/converging  nozzle. 

Boundary  Layer  Generation 

The  throat  area  and  the  nozzle  exit  area  have  been  calculated  using  a  mass  flow  rate 
that  incorporates  power  losses  by  the  engine.  While  these  results  are  suitable  for  designing  an 
engine,  it  is  important  to  calculate  the  boundary  layer  associated  with  the  flow,  even  though  the 
losses  associated  with  boundary  layer  development  are  typically  only  1-10%  of  total  flow  area. 
For  larger  systems,  this  is  significant  because  only  approximately  90%  of  the  flow  remain 
undisturbed  as  it  is  accelerated  and  expanded  through  the  engine.  The  boundary  layer  is 
calculated  by: 

5b  =  0.16x/(Rex)1/4.  (17) 

For  turbulent  flow  the  Reynolds  number  (Re,)  is  given  by: 

Rex  =  vx/v  where  vx=  [(2gcyRT1)/(Y+ 1)]iy2.  (18) 

The  thermal  expansion  correction  factor  is  given  by: 

8,h  =  a(AT)L.  (19) 

In  Equation  15  alpha  (a)  is  material  specific  and  AT  is  experimentally  obtained.  These  two 
parameters  are  then  multiplied  by  nozzle  length  (L).  Nozzle  length  for  STP  would  be  the  length 
along  the  thrusting  axis  starting  at  the  absorber  exit  plane  ending  at  the  engine  exit  plane.  Note 
the  length  of  the  nozzle  must  be  truncated  to  discourage  flow  separation  and  shock  formation. 

The  final  correction  factor  relates  the  change  in  the  nozzle’s  effective  circumference  due  to 
pressure  effects.  This  last  correction  factor  also  includes  stress  effects: 


c  =  (P2-P3)r/t.  (20) 

In  this  equation,  sigma  (a),  is  a  function  of  throat  and  exit  pressure  in  conjunction  with  thickness 
and  nozzle  throat  radius.  The  radius,  r,  is  the  actual  radius  without  any  correction  factor.  The 
final  radius  will  be  calculated  from  the  sum  of  these  losses.  The  thickness  of  the  material,  t,  is 
also  a  known  parameter.  It  is  determined  by  the  availability  of  the  material  used  (typically 
tungsten,  rhenium,  molybdenum)  and  material  requirements  (see  Equation  15).  The  value  that  is 
calculated  for  a  can  then  be  related  to  an  actual  lessening  of  circumference  in  Equation  21. 


8cir  ”  o  L  /  E 


(21) 


where  E  is  the  modulus  and  is  a  material  specific  value.  The  change  in  the  circumference  can 
then  be  converted  to  a  change  in  diameter.  The  total  correction  in  the  throat  diameter  is  found  by 
the  summation  of  all  of  the  correction  values: 


8tot  =  5cir  "*■  5b  =  6th-  (22) 

Equation  22  may  now  be  applied  to  the  throat  diameter.  This  correction  equation  includes  the 
boundary  layer,  thermal  expansion  and  pressure  influence. 

While  it  is  not  much  of  technical  challenge  to  use  modified  equations  to  produce  an  engine 
that  meets  specific  mission  requirements,  actually  fabricating  and  testing  the  designed  engine 
may  prove  to  be  difficult.  For  low-level  thrust  requirements  typically  the  “nozzle”  employed  has 
been  a  bored-out  hole  aft  of  the  absorber.  This  method  works  reasonably  for  chamber  tests 
(order  of  1  x  1C6  torr),  but  is  not  optimal  for  space-rated  hardware.  Special  fabricating 
procedures  must  be  practiced.  Precise  equipment,  preferably  machinery  used  in  conjunction  with 
a  CAD  package,  is  mandatory  to  meet  the  tolerances  required  for  a  STP  system. 

Once  the  engine  has  been  built  to  specifications,  testing  may  be  initiated.  In  the  ideal 
scenario,  thrust  measurements  taken  during  the  tests  will  reinforce  numerical  estimates. 
Unfortunately,  there  is  no  foolproof  method  for  achieving  the  same  numbers  experimentally  as 
analytically.  It  is  evident  that  while  these  equations  help  reduce  engine  performance  uncertainty, 
they  promise  no  guarantee  for  performance  results.  The  most  equitable  method  for  ensuring 
rocket  performance  is  to  test  with  the  correction  factor  iterations  until  discrepancies  can  be 
identified  and  corrected. 

Each  engine  that  will  be  used  for  STP  will  change  with  changed  mission  requirements.  For 
example,  where  one  mission  may  require  a  simple  orbit  transfer  of  a  satellite,  another  may  be 
specified  to  execute  a  deep  space  probe.  Chemical  propulsion  systems  are  currently  reaching 
mission  limitations  due  to  high  fuel  requirements  and  propellant  storage  limitations.  If  deep 
space  missions  are  to  continue,  an  earth-to-orbit  chemical  propulsion  system  with  a  high  Isp 
upper  stage  must  be  used  to  lower  launch  costs  and  reduce  lift-off  weight.  The  solar  group  plans 
to  test  a  variety  of  engines  to  meet  USAF  propulsive  upper-stage  needs  and  share  the 
technology  they  develop  with  the  commercial  launch  vehicle  industry. 

Conclusions 

Most  of  the  engines  tested  in  the  solar  laboratory  have  been  designed  and  fabricated  by 
government  contractors  as  deliverables  in  SBIR  (small  business  innovative  research)  contracts. 
After  an  engine  is  produced,  it  is  shipped  to  the  Solar  Lab  where  it  undergoes  on-sun  testing  with 
a  ground-based  concentrator.  A  data  acquisition  system  monitors  characteristic  pressures  and 
temperatures  and  relates  them  to  performance  values.  (For  typical  test  results  see  Figure  4.) 


Figure  4.  Relationship  between  Thrust,  Temperature  and  Isp 

Some  members  of  the  solar  group  also  design  and  build  engine  components  and  systems 
fortesting.  The  work  presented  in  this  paper  will  be  used  by  the  members  in  the  solar  group,  and 
will  be  compared  with  contractor  engine  studies.  While  most  of  the  engine  work  done  to  date  has 
been  on  the  blackbody  absorber  type  of  engine,  the  equations  presented  in  this  paper  can  also 
be  applied  to  other  engine  types  (Reticulated  Vitreous  Carbon  Thruster  (RVCT),  Beam- 
Fractionating,  et  al.). 

Tests  conducted  in  the  near  future  will  validate  AFRL’s  ability  to  improve  current  engine 
configurations.  These  tests  will  demonstrate  the  STP’s  ability  to  operate  for  extended  periods  of 
time  under  high  temperatures  and  to  continually  meet/exceed  IHPRPT  goals.  IHPRPT, 
Integrated  High-Payoff  Rocket  Propulsion  Technology,  is  the  means  by  which  in-house  efforts 
receive  funding.  Performance  goals  must  be  met  to  continue  IHPRPT  funding.  By  continually 
improving  the  design,  engines  are  continuously  and  consistently  out  performing  their 
predecessors. 
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